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pL( ' Abstract. Galaxy-galaxy lensing (GGL) measures the 2-point cross-correlation between galaxies and mass in the 

, Universe. In this work we seek to generalise this effect by considering the third-order correlations between galaxies 

^s^j ■ and mass: galaxy-galaxy-galaxy lensing. Third-order correlations in the cosmic shear field have recently been 

reported in the VIRMOS-DESCART and CTIO surveys. Such data should also be ideal for measuring galaxy- 
y—i . galaxy-galaxy lensing. Indeed, the effects of these higher-order correlations may have already been detected in 

^ ' recent studies of galaxy-galaxy lensing. Higher-order cross-correlation functions contain invaluable information 

about the relationship between galaxies and their mass environments that GGL studies alone cannot detect, 
in . In this paper we lay out the basic relations for third-order cross correlations and their projections and introduce 

a new set of scale dependent third-order bias parameters. We define three new observables: two galaxy-shear-shear 
correlation functions, G±, and a galaxy-galaxy-shear correlation, Q. We relate these to the various projected cross- 
bispectra and give practical estimators for their measurement. We note that the observational signature of these 
correlators is simply the excess shear-shear correlation measured about foreground galaxies (for G±) and the 
average tangential shear around foreground galaxy pairs (for Q). These quantities are no more than second order 
' in the shear and so should be more easily measurable than the shear 3-point correlation. Finally we derive 

expressions for the third order aperture mass statistics in terms of both the cross-bispectra and the real-space 
correlation functions. Such statistics provide a very localized measurement of the bispectra, thus encapsulating 
essentially all of the available third-order information, while remaining easily obtainable from observations of 
3-point cross-correlation functions. In addition we find that utilising aperture statistics has the further benefit 
that they measure only the connected part of the third order correlation. 
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1. Introduction 

Gravitational lensing offers the possibility of studying the statistical properties of the mass distribution in the Universe 
without referring to the relation between the mass and the luminous tracers, like galaxies. The weak lensing effect (see 
Mellier 1999; Bartelmann & Schneider 2001, hereafter BS01; Refregier 2003; van Waerbeke & Mellier 2003 for recent 
reviews) of the large-scale structure, called cosmic shear (see, e.g., Blandford et al. 1991; Miralda-Escude 1991; Kaiser 
1992; Jain & Seljak 1997), has been detected by a number of groups (e.g. van Waerbeke et al. 2000, 2001, 2002; Bacon 
et al. 2000; Kaiser et al. 2000; Wittman et al. 2000; Maoli et al. 2001; Refregier et al. 2002; Hammerle et al. 2002; 
Brown et al. 2003; Jarvis et al. 2003). Most of these cosmic shear measurements have concentrated on second-order 
shear statistics, but quite recently, third-order cosmic shear measurements have been reported as well (Bernardeau et 
al. 2002; Pen et al. 2003; Jarvis et al. 2004). 

Weak lensing can also be used for measuring the relation between mass and galaxies. The measurement of this 
galaxy-galaxy lensing (GGL) was first attempted by Tyson et al. (1984), and its first detection was published by 
Brainerd et al. (1996). Since then, quite a number of GGL measurements have been reported, in the HST Medium 
Deep Survey (Griffiths et al. 1996), the Hubble Deep Field (DelL Antonio & Tyson 1998; Hudson et al. 1998), and 
more recently, in the RCS survey by Hoekstra et al. (2001, 2002), the COMBO-17 survey (Kleinheinrich et al. 2004) 
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and the Sloan Digital Sky Survey (Fischer et al. 2000; McKay et al. 2001; Guzik & Scljak 2002; Sheldon et al. 2004, 
Seljak et al. 2004). GGL measures the two-point correlation function (2PCF) between galaxy positions and shear. 

On small scales, GGL measures the lensing effect of the dark halo in which the luminous galaxy is embedded. It 
therefore provides a probe of the density profile of galaxy halos which can be compared to the predictions of the CDM 
model. In particular, the properties of dark matter halos, such as virial mass and concentration, can be studied as 
a function of galaxy type (early vs. late), environment (e.g., have the halos of galaxies in groups and clusters been 
stripped by the tidal interactions), and redshift. On scales larger than ~ 200/i~ 1 kpc, the GGL signal is no longer 
dominated by the halo of galaxies, but receives a substantial contribution from the matter in which galaxies are 
embedded, namely groups and clusters. From this signal, one learns about the environments of galaxies. For example, 
the density-morphology relation can be traced back directly to the underlying mass density, instead to the galaxy 
number density. 

Schneider (1998) pointed out that the correlation between galaxy positions and shear provides a direct measure of 
the bias factor of galaxies; in particular, the scale dependence of the bias factor can be probed directly (van Waerbeke 
1998). Indeed, on scales ^ 2/i~ 1 Mpc, GGL measures a combination of the galaxy-mass correlation coefficient r and 
the bias factor b. These two functions can be obtained individually when the GGL signal is combined with the cosmic 
shear signal, as shown by Hoekstra et al. (2002). In fact, Hockstra et al. (2002) have measured r and & as a function 
of scale. Future measurements of these quantities as a function of scale, redshift and galaxy type, using the upcoming 
wide-field imaging surveys, will provide invaluable input for the interpretation of galaxy redshift samples. In addition, b 
and r are functions which depend on the formation and evolution of galaxies; they can therefore be used as constraints 
for those models. 

As already mentioned, a non- vanishing cosmic shear three-point correlation function (3PCF) has been detected, and 
future surveys will measure it with great precision. These surveys will certainly also measure higher-order correlations 
between galaxies and matter, generalizing the GGL measurements. Although the detailed physical interpretation of 
such higher-order cross correlations remains to be investigated, it is obvious that they contain valuable information 
concerning the relationship between galaxies and their mass environments beyond that which is measurable in GGL. 
In the language of the halo occupation distribution (Bcrlind & Weinberg 2002), one can think of these higher-order 
functions as probing the moments of an occupation probability P{N\M) that a halo of virial mass M contains N 
galaxies of a particular type. Such information is of crucial importance since it places powerful constraints upon 
models of galaxy formation. 

We save the detailed interpretation of higher-order cross correlations for a forthcoming paper, where we explore 
such themes in the context of the halo model. In this work we lay out the basic relations for the third-order galaxy-mass 
correlations. After introducing the projected mass and galaxy number densities in Sect. 2, we define the third-order 
bias factor and two third-order correlation coefficients in Sect. 3. The relation between the projected bispectra and 
the three-dimensional bispectra is given in Sect. 4. Practical estimators for the galaxy-galaxy-shear and the galaxy- 
shear-shear correlation functions are provided in Sect. 5. We will argue that these third-order galaxy-mass correlations 
have already been measured and published. In Sect. 6, the relation between the correlation functions and the cross- 
bispectra are derived, and a consistency relation for the two galaxy-shear-shear 3PCFs is obtained that should be 
satisfied provided the shear is a pure E-mode field. In Sect. 7 we define aperture measures of the projected matter 
and galaxy densities and consider their third-order statistics; in particular, on the one hand they are related to the 
projected bispectra, on the other hand, they can be directly calculated from the respective correlation functions. We 
discuss our results in Sect. 8. 

Menard et al. (2003) have considered a different approach to measure third-order correlations between the mass and 
the galaxy distribution, employing the magnification bias of background galaxies. Generalizing an effect that has been 
observed at the two-point statistical level, where correlations between high-redshift QSOs and low redshift galaxies are 
observed and which are interpreted as being due to the magnification caused by the large-scale matter distribution of 
which the foreground galaxies are biased tracers (see Dolag & Bartelmann, 1997; Bartelmann & Schneider 2001, and 
references therein), Menard et al. (2003) considered the excess of distant QSOs around pairs of foreground galaxies. 
For linear deterministic bias of the galaxy distribution, the ratio of this third-order statistics and the square of the 
QSO-galaxy correlation function becomes essentially independent of the bias, as well as of the shape and amplitude of 
the power spectrum of density fluctuations. This result is very similar to the corresponding ratio of third-order cosmic 
shear statistics and the square of the shear dispersion (Bernardeau et al. 1997; Schneider et al. 1998; van Waerbeke et 
al. 1999). 

2. Projection of density and galaxies 

The gravitational lensing effect of the inhomogeneous matter distribution in the Universe is described in terms of an 
equivalent surface mass density k(0), which is obtained by projecting the three-dimensional density contrast, 8 7 of the 



Peter Schneider and Peter Watts: Galaxy-galaxy-galaxy lensing 



3 



matter along the line-of-sight (see BS01, and references therein). If we consider sources at comoving distance w, this 
surface mass density is given by 

2c 2 J Jk(w) a(w') 

where Jk{w) is the comoving angular-diameter distance corresponding to comoving distance w, /if(u>) — w for flat 
Universes, 6 = Ap/p is the relative density contrast of matter, Hq is the Hubble constant, c the velocity of light, 
a(w) = 1/(1 + z) is the cosmic scale factor, normalized to unity today, and f2 m is the density parameter in matter. For 
a redshift distribution of ('background') sources with probability density p z (z)dz — p w (w) dw, the effective surface 
mass density becomes 



(2) 



<f>\ (a < \ (a \ 3g o 2 "m P h , ( w , , S(f K (w)0,w) 
k{0) = J dw p w {w) k{0,w) = ^ c 2 J &wg{w)f K (w) , 

with 

g{w) = / dw p w (w ) — , , (3) 

Jw Jk{w ) 

which is the source-redshift weighted lens efficiency factor D^ s /D s for a density fluctuation at distance w. The quan- 
tity Wh is the comoving horizon distance, obtained from w(a) by letting a — > 0. The shear, i.e. the projected tidal 
gravitational field which can be measured as the expectation value of image ellipticities of the galaxy population, is 
obtained from the projected surface mass density k in the usual way: one defines the deflection potential ip(0) which 
satisfies the two-dimensional Poisson equation V 2 V> = 2k, in terms of which the Cartesian components of the shear 
read 71 = (ip,u — ^,22) /2, 72 = tp, 12, where indices separated by a comma denote partial derivatives. 

In analogy to the matter density contrast S, we define the number density contrast S g (x, w) of galaxies as 

n(x, w) - fi(w) 

5 K (x,w):=-± — _ 4 
n(w) 

where n(x, w) is the number density of galaxies at comoving position x and comoving distance w (the latter providing 
a parameterization of cosmic time), and n(w) is the mean number density of galaxies at that epoch. Since the galaxy 
distribution is discrete, the true number density is simply a sum of delta-functions. What is meant by n is that the 
probability of finding a galaxy in the volume dl^ situated at position x is n(x) dV. We consider now a population 
of ('foreground') galaxies with spatial number density n(x,w). The number density of these galaxies on the sky at 
6 is then N(0) = J dw v(w)n(fk(w)0,w), where v(w) is the redshift-dependent selection function, describing what 
fraction of the galaxies at comoving distance w is included in the sample. The selection function v(w) can differ 
between different types of foreground galaxies and thus depends on the sample selection. The mean number density 
of these galaxies on the sky is N — J dw v(w) n(w); the redshift distribution of these galaxies, or more precisely, their 
distribution in comoving distance therefore is pf(w) — v(w)n(w)/N, thus relating the selection function v{w) to the 
redshift distribution. Using the definition one then finds for the number density of galaxies N(0) on the sky and 
their fractional number density contrast k s (0), 

N(0) =N 1+ f dw Pi (w)6 g (f K (w)0,w) ; K g (0) := [N(0) - N] /N = J dw p { (w) 8 & {f K {w)0,w) . (5) 



3. Definition of bias and correlation coefficients 

3.1. Second-order statistics 

In its simplest form, the relation between the density contrast 5 and the fractional number density contrast 6 g of 
galaxies is described by the bias factor b, in the form S g — bS. This picture of linear deterministic biasing is most likely 
too simple a description for the galaxy distribution, at least on small spatial scales where the density field is non-linear. 
Instead, one defines the bias parameter in terms of the power spectra of the galaxy and matter distribution. Let 5(k) 
and 5 g (k) be the Fourier transforms of the density and galaxy contrast, respectively; their power spectra P$$ and P gg 
are defined by the correlators 

(S(k,w)5*(k',w)\ = (2TT) 3 S D {k-k')P ss {\k\,w) ; (& g (k,w)%(k',w)\ = {2jr) 3 5 D {k - k') P gg {\k\,w) , (6) 

where <5d denotes Dirac's delta-'function', and we have allowed for an evolution of the power spectra with redshift 
or comoving distance w. The occurrence of the delta function is due to the assumed statistical homogeneity of the 
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random matter and galaxy fields, and the fact that the power spectra depends only on the modulus of k is due to 
their statistical isotropy. The bias parameter b(k,w) is defined by the ratio of these two power spectra, 

P gg (k,w) =b 2 (k,w)P ss (k,w) , (7) 

which agrees with the previous definition in the case of linear deterministic biasing, but is far more general. In 
particular, we allow for a scale- and redshift dependence of b. Next, we define the cross-power spectrum P$ g (k,w) 
through 

(s(k,w)5*(k',w)} = (2ir) 3 6 B (k-k')P Sg (\klw) ;. (8) 
The cross-power spectrum is related to the power spectrum of the matter density by 

Ps g (k,w) = b(k,w)r(k,w) Pss(k,w) , (9) 

where we have defined the galaxy-mass correlation coefficient r(k, w), which can also depend on scale and redshift. In 
the case of linear deterministic biasing, r = 1. 

3.2. Third-order statistics 

We now generalize the foregoing definitions to third-order statistics. The bispectra of the matter and galaxy distribu- 
tions are defined through the triple correlators 

l(fci) J(fe 2 ) <5(fc 3 )) = (2 7 r) 3 fc(fc 1 +k 2 + k 3 )B sss (k 1 ,k 2 ,k 3 ;w) , 
4(fci) S s (k 2 ) S g (k 3 )) = (2tt) 3 fo(fc a + fc 2 + fe 3 ) B ggs (k 1 ,k 2 , k 3 ; w) , (10) 



where the delta function ensures that these triple correlators vanish unless the three fc-vectors form a closed triangle 
- this property again follows from the statistical homogeneity of the density fields. Furthermore, statistical isotropy 
causes the bispectrum to depend only on the length of two of its fc-vectors and the angle they enclose. Finally, parity 
invariance of the density fields requires that the bispectra are even functions of this angle. In a similar spirit to Eq. JHJ 
we can define the cross-bispectra 

<5(fci)<5(fc 2 )4( fe 3)) = (27r) 3 S (ki + k2 + k 3 )Bss g (kx,k2;k 3 ;w) , 

S s (ki) 8 g {k 2 ) S(k 3 )) = {2nY <5 D (fci +k 2 + k 3 ) B ggS {k u k 2 : k 3 ; w) . (11) 



where we have introduced the notation (ki, k 2 ; k 3 ) to indicate the symmetry with respect to interchanging the first 
two arguments, meaning, e.g., that B gg s(ki, k 2 ; —k\ — k 2 ) — B gg s(k 2 , kx; — k\ — k 2 ). 

In the case of linear deterministic biasing, the various bispectra would be related through simple relations, e.g., 
Bs gg = b 2 Bsss- However, in the more general (and realistic) case, these relations are more complex. We define the 
third-order bias parameter b 3 and the two galaxy-mass correlation coefficients r\ and r 2 through the relations 

Bg g g(fci,fc 2 ,fc3; w) = bl(kx,k 2 ,k 3 ;w) Bsss{ki,k 2 ,k 3 ;w) , 

B gg s{ki,k 2 ]k 3l w) = bl(k 1 ,k 2} k 3 ;w)r 2 (k 1 ,k 2 ;k 3 ;w) B S ss(ki 7 k 2 ,k 3 ;w) , (12) 
Bs6 s {ki,k 2 ;k 3 ] w) = b 3 (k 1 ,k 2 ,k 3 ;w)r 1 (k 1 ,k 2 ;k 3 ;w) B S ss{ki,k 2 ,k 3 ;w) . (13) 

The bias and correlation coefficients satisfy the same symmetries as those of the bispectra. In the case of linear 
deterministic biasing, b 3 = b and r\ = 1 = r 2 . These higher-order biasing and correlation parameters depend on the 
formation and evolution of galaxies and their relation to the underlying density field; in general, they encapsulate 
information concerning the non-Gaussian nature of the bias. As we show in the next section, the correlation of lensing 
shear with the positions of galaxies on the sky provides a direct means for measuring (integrals of) these higher-order 
functions. 



4. Projected power- and bispectra 

The lensing effect and the position of galaxies on the sky are both described by projections of the density field S 
and the galaxy field 5 g along the line-of-sight, as shown in Sect.[5] The projected fields share the properties of the 
three-dimensional distributions as being homogeneous and isotropic random fields; hence, their second- and third- 
order properties can as well be described by power and bispectra. The relation between the power spectra of projected 
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quantities to that of the three-dimensional distribution is given by Limber's equation in Fourier space, as has been 
derived by Kaiser (1992): Let 

Ki(0) = J dw qi (w) St (f K (w)0, w) (14) 

be projections of the 3-D fields Si, where the qi{w) are weight functions depending on the comoving distance. Provided 
these weight functions do not vary appreciably over comoving scales on which the power spectrum is significantly 
different from zero, the (cross) power spectrum reads 

qi(w)q 2 (w) ( I 



Pi2(£)= dw^f^Psl-— w . (15) 

Hence, the 2-D power at angular scale 2ir/£ is obtained from the 3-D power at length scale Jk{w) (2tt/£), integrated 
over w. For the power spectra of the projected mass density k and the galaxies K g , and the cross-power of these two 
quantities, we therefore find, by comparing Eqs. Ijl4(l and l|15l) with (J2J) and JHJ), 

P KK {£) = [ Wh dw^-P 5S (k;w) ; P^i) = f dw ifM b 2 (k; w) P 5S (k; w) ; 

4c Jo a \ w ) J Jk( w ) 

PM = / ^ f^P^\ b (fc; w) r (fc; w) Pss (fc; w) , (16) 

where k = £/ fx(w) and where we have used the definition of the bias and galaxy-mass parameters given in Sect. 13. ll 
The quantity P KK (£) is the power spectrum probed by second-order cosmic shear measurements whereas P sg (£) is the 
power spectrum of the galaxy distribution [with the selection function specified by pi(w)] on the sky, i.e., the Fourier 
transform of the angular 2PCF of galaxies. The cross-power P Kg (£) is the quantity that is measured in galaxy-galaxy 
lensing observations; see Hoekstra et al. (2002). 

The bispectrum of the projected quantities (|14fl is defined in terms of the triple correlator 

(ki(£i)« 2 (€2)k 3 (^)) = (2vr) 2 ,5D^i+^+^)foi23^i^2,4) , (17) 
where the relation between the bispectrum of the projected quantities k and the 3-D bispectrum is given by 

t>123 (*1,«2,«3) = dw - r — #123 7 , v t 7 v t , ^ W , (18) 

J Ik\ w ) \Jk{w) f K (w) f K [w) J 

which is valid under the same conditions required for the validity of (|15fl . i.e., the weight functions qi(w) should not 
vary appreciably over scales on which the bispectrum is significantly non-zero. Inserting into this relation the weight 
functions q appropriate for the surface mass density k and the projected number density K g , one obtains the bispectra 



b K K,K,{£l, £2: £3) 


8c b J /Ar(w)a d (w) 






(19) 


b KKS (£ u £2;£ 3 ) 


= " 4 / dw J \ >; 2 ) > [b 3 r 1 B S ss](k 1 ,k 2 ;k 3 ;w) =: 
4c 4 J f^(w)a z (w) 


[b 3 r 1 b K ^ K ] (£1, 


£2-^3) , 


(20) 


b EEK ,{£i,£2',£z) 


3H^Q m f g{w)p 2 (w) r 2 -, 
= n 2 / dw ' [b 3 r 2 B sss \(k 1 ,k2;k 3 ;w) =: 

2c J JK{ w ) a \ w ) 


[blr 2 b KKK ] (£1,. 


£2^3) 


(21) 


b ges (£i,£ 2 ,£3) 


= fdw §^ [blB 5SS ](k ll k 2 ,k 3 ;w) =: [b 3 3 b KKK ] (£ h £ 2 , 
J Jk\ w ) 


£3) , 




(22) 



where again the vectors k L = £ij Jk(w), and we have defined the projected bias factor 63 and correlation coefficients fi, 
f 2 , which depend, in addition to the angular wave vectors £i, on the redshift distribution of foreground and background 
galaxies. The correlations are strongest if the redshift distribution of the foreground galaxies matches the most efficient 
lensing redshift for a given source redshift distribution - see Hoekstra et al. (2002) for a similar discussion in the context 
of second-order galaxy-mass correlations. 



5. Estimators 

The basic observables for the second-order statistics obtained from observational data are the 2PCFs. Hence, the 
angular 2PCF of galaxies w(9) = (K g (<p)K g (<p + 9)) is the Fourier transform of P gg (£). The second-order statistics 
of the projected matter density are probed by the shear 2PCFs £±(0). They are defined in the following way: let 
7 C = 71 + i72 denote the complex shear in Cartesian coordinates. Given a direction ip, we define the rotated shear 
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Fig. 1. Geometry of the galaxy-shear-shear correlation, G± (#1, 1)2, 03 ) (left panel), and the galaxy-galaxy-shear cor- 
relation, Q §2i 03 ) (right panel). Note that the sign of the angle <fo = ip2 — ifi is important. 

~f(8',(p) = — e _2lv 7 c (0). This definition implies that 7(6; ^1) = j{6; P2) e 21 ^ 2 ^ 1 ^. The real and imaginary part of 
7(0; <p) are the tangential and cross component of the shear relative to the direction tp (see, e.g., Crittenden et al. 
2002; Schneider et al. 2002). Hence the shear 2PCFs are written as 



where the average is over all pairs of points with angular separation (9, and tp is the polar angle of their separation 
vector 9. The imaginary parts of £± vanish due to parity invariance (e.g., Schneider 2003). For a set of observed 
galaxies, the shear in (|23|l is replaced by the image ellipticities, which are an unbiased estimator of the shear (we 
neglect here the difference between shear and reduced shear; see BS01). The shear 2PCFs are related in a simple way 
to the power spectrum P KK (£) - see Kaiser (1992). 

The cross-power is probed by the galaxy-galaxy lensing signal, which correlates the tangential component of the 
shear at the location of a background galaxy with the position of a foreground galaxy. We define this correlation 
function as 



(e.g., Hoekstra et al. 2002), where ip is, as before, the polar angle of the connection vector 9, and the imaginary part 
of Q24JI vanishes due to parity invariance. A practical estimator for (jt) is obtained by averaging the tangential shear 
component of the background galaxies in all foreground-background pairs with separation 8. 

Similarly, the third-order statistics are best probed from observing the 3PCFs. The angular 3PCF of galaxies is one 
of the standard ways to characterize the non-Gaussian properties of the galaxy distribution on the sky, and standard 
methods for estimating it are known (see, e.g., Peebles 1980). The bispectrum b KKK is probed by the shear 3PCF, 
as defined in Schneider & Lombardi (2003; see also Takada & Jain 2003a; Zaldarriaga & Scoccimarro 2003), and its 
relation to the bispectrum has been derived in Schneider et al. (2004). We next derive expressions for the correlation 
functions related to the cross-bispectra. 

5.1. Galaxy-shear-shear correlation 

Consider a triplet of points, where two background galaxies are located at 9\ and 92, and a foreground galaxy is 
located at #3 (see Fig. 1, left). We define the following galaxy-shear-shear 3PCFs 



£ f (0) = <7to¥>)7*(0 + 0;p)> ; £-(*) = (7(6 *>)7(* + 0;v)> > 



(23) 



(7t> W = </%(<£) 7(0 + e-,v)) 



(24) 



G+(tfi,i? 2 ) = G+(0i,0 2> 03) = < 7 (0i; </?i)7*(0 2 ; <Pi) «g(0 3 )> 
G_(0i,0 a ) = G_(0i,i^,&) = <7(0i;¥>i)7(02;¥>2)/s g (03)) , 



(25) 
(26) 



where the separation vectors between the background galaxies and the foreground galaxy are #i = 9i — 9% for 
i = 1,2, and ipi is the polar direction of the separation vector These correlation functions depend on the moduli 
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of the separation vectors 0i and the angle 03 enclosed by them. In order to find a practical estimator for these 
cross-correlations, it is useful to consider the modified correlators 

G + (#i,# 2 ,(l )3 ) = ±( 1 (e 1 ;<p 1 )~f*(e 2 ;<p 2 )N(d s )) , G_(tfi,tf 2 , 3 ) = 1 (7(61; <Pi) 7(^2; ^2) JV(6 3 )> . (27) 
Using A^/TV = 1 + Kg, these modified correlators become 

G+(0i,02,0 3 ) = G + (0i,0 2) ^3) + (7(6'i;^i)7*(^;^2)) = G+(0i,02,03)+ 1 e + (A&)e 2i ^ , (28) 

G_ (0i,0 2 ,03) = G-(0i,02,^3) + (7(ei;yi)7(g2;y 2 )) = g-(0i,02,^3) + e-(Ag) ^ 2e ' 03/2 (A J ) 1 4 6 ' 03/2 ^ , (29) 

where AO = |0 2 — i?i | is the separation between the two background galaxies. The phase factors with which the shear 
2PCFs are multiplied are obtained by rotating the shears from the direction pointing towards the foreground galaxy 
to that of the separation vector of the two background galaxies. If the polar angle of the latter is denoted by ifA, then 
7(0;; (pi) = j(6i; <^a) e 2l< ' VA ~ Vi \ from which the expressions for the phase factors follow. Hence, the modified galaxy- 
shear-shear correlators are given by the sum of the (reduced) galaxy-shear-shear correlators G± (which is proportional 
to the cross-bispectrum b KKg ) plus the 2PCF of cosmic shear (|23[1 . modified by a phase factor that accounts for the 
fact that the projection directions of the shear are ifi, and thus do not correspond to the separation vector between 
the background galaxy pair. 

We can obtain practical estimators for these cross 3PCFs, as follows: first, one defines bins in (0i, 2 , 3 )-space. 
Then, for each triplet of two background galaxies and one foreground galaxy whose separation falls inside a given bin, 
one adds the product of shears, so that 

G+(0 1) 02,0 3 ) - — — 7(015 ^i)7*(^; V2) ; G_(0i,0 2 ,0 3 ) - 7(011 ViMhit*) • (30) 

Nt "V let triplet Nt "P lct triplet 

These estimators can be easily calculated from the catalog of galaxy positions and ellipticities, e.g., using the method 
of Jarvis et al. (2004) for finding triplets. From the estimator above, the contribution of the shear 2PCF according to 
l|29|) has to be subtracted in order to obtain the reduced galaxy-shear-shear 3PCF G±. Note that observationally the 
galaxy-shear-shear correlation is simply the excess of the shear 2PCF measured around foreground galaxies. 

5.2. Galaxy-galaxy-shear correlation 

Next we consider the 3PCF between the shear of one background galaxy and the positions of two foreground galaxies. 
Let the former be at position #3, the latter at positions 0\ and 62, then we define the reduced galaxy-galaxy-shear 
correlation function as 

5(01, 02) = 5(01, 02, 03) = /«g(0l)«g(02)7 (03; ^1±^\ \ , (31) 

where, as before, 6i = 9 3 + 0^, and ipi is the polar angle of the vector 0j, for i — 1,2 (see Fig. 1, right). The shear 
in l|31|) is projected along the line which bisects the angle between the two foreground galaxies, as measured from the 
background galaxy. As before, for obtaining a practical estimator it is useful to first define a modified 3PCF, 

<?(0i, 2 ) = 5(0i, 2 , 3 ) = ij (^N(6 1 )N(e 2 ) 1 (d 3 ; fl±^lj J . (32) 

Using again N/N = 1 + k s , and applying the transformation law for rotating shears, this can be written as 

5(0i, 2 ) - 5(01, 02, 3 ) = 5(01, 02, 3 ) + <7t> (0i) + ( 7 t) (0 2 ) e i<fe , (33) 

where the last two terms are the two-point galaxy-shear correlation functions l|24|l . obtained from rotating the shear 
in the direction of the two foreground galaxies. Hence, the correlator Q is the sum of the reduced galaxy-galaxy-shear 
correlation function plus the two GGL contributions which would be present even for the case of purely Gaussian 
density fields. A practical estimator for Q is obtained by finding triplets of galaxies that fall in a given bin, and then 
summing over the shears of the background galaxies, 

5 = ^V 7 ^;^ ± ^V (34) 

triplet tr tTet V 2 / 

The signature of the galaxy-galaxy-shear correlation is therefore just the excess shear measured about pairs of fore- 
ground galaxies. 
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Note that integral measures of these three-point cross-correlation functions have probably been measured already. 
McKay et al. (2002) have demonstrated in their measurement of GGL from the SDSS that the GGL signal is stronger, 
and extends to much larger separations, for foreground galaxies that are located in regions of high galaxy density. 
This detection provides a correlation between the GGL signal and the number density of galaxies, and therefore an 
integral over the 3PCF Q. Furthermore, the galaxy-shear-shear correlation seems to be in the cosmic shear analysis 
of the COMBO-17 fields by Brown et al. (2003), where they find a stronger-than-average cosmic shear signal in the 
A901 field, and a weaker-than-average cosmic shear signal in the CDFS, which is a field selected because it is rather 
poor in brighter galaxies. 

The measurement of these galaxy-mass 3PCFs is expected to be considerably easier than measuring the shear 3PCF 
itself. Note that Q is first order, and G± are second order in the shear. For the same reason that GGL measurements 
are easier to obtain than second-order cosmic shear measurements, the former being just first order in the shear, one 
therefore expects that Q and G± are more straightforward to measure from a given data set. fn particular, every wide- 
field survey in which the shear 3PCF can be measured should yield a significant measurement of these cross-correlation 
functions. 



6. Relation to the bispectrum 

It is relatively straightforward to write the correlation functions G± and Q in terms of the (projected) bispectra defined 
in Sect. 01 The complex shear is related to the convergence in Fourier space so that %(£) — k(£)e 2l P, where (3 is the 
polar angle of the vector t. Using this, along with the rule for rotating shears, one can re-write equation l|25|) as 



(35) 



Similar relations hold for the correlation functions defined by Eqs. Q2(j[l and (|31|l . Next we split the angles tp i and $ 
into their mean and their difference, by writing y>\ = Q — cj)z/2, <p 2 — C + 03/2 and f3\ = rf — tp/2, (3 2 = if + "0/2, 
where ip represents the angle contained by £\ and £ 2l and 03 is the same as in previous sections. Then, we write 
i\ ■ i?i + £ 2 ■ i?a = A cos(?7 — v), where rj = rf — Q, and 

A 2 = i\d\ + £\d\ - 2li£ 2 did 2 cos(0 3 - , e 2iu = [2Wn?2 + (4#i) 2 e 1 ^ 3 ^ + {l 2 ti 2 f e^-V*)] . ( 36 ) 
Splitting up the integrals in l|35|) into polar coordinates in £-space, the ^-integral can be performed, yielding 

G + (^ 2 M=J ^/f^ e- 2i( ^ MA) . (37) 

In a similar manner, we find for G_ 

G^ u # 2 ,fo) = J ^2 J ^e^-^+^^'e^^^e 2 ^-^)^^,^;-^-^) (38) 

dttii fd£ 2 e 2 r dt/j ilv 

WlWlW e MA) , (39) 

and for Q 



with \£\ 2 = £ 2 + £\ — 2l\£ 2 cosip. The above results link directly the cross-correlation functions in real space to their 
equivalent projected bispectra. In a similar fashion, one may also derive a set of inverse relations, i.e expressing the 
bispectra in terms of Fourier integrals of the correlation functions, beginning with the definition of the bispectrum in 
Eq. (|f 7|l . The resulting expressions, which are derived in the same way as above, are 

b KKS (£x,£ 2 ;-£x-£2) = J d 2 ^J d^e^-^+^-^Je-^-^e-^-^G-^,^) (41) 
= J d 2 tfr J d 2 d 2 e K£i-0i+£2-&2) e -HPi- V i) e 2i(A- V2 ) G+iWi) . (42) 

(43) 
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Using the same transformations of the polar angles as before, one more integration can be carried out, resulting in 
W*i>*a,V) = ^ J d ^i #1 J d #2 ^2 Jdfo e 4i " h{A) G_(0i,02, 03 ) (44) 
= 2tt /dt?n?i /"d0 2 2 /d0 3 e 2i ^-«Jo(A)G + (^,i?2,</»3) ■ (45) 



In a similar way, we can express the cross-bispectrum b ggK in terms of Q, yielding 

r r r (L^/ 2 + i->e-^' 2S \ 2 

b gSK {£i,£ 2 ^) = 2nJ &0i t?i J d0 2 2 J d0 3 ^ ^ e- 2iu J 2 (A) Qtfu^fa) . (46) 

Using l|35|> and replacing the cross-bispectrum in favour of G_, using l|41|) . one can express G+ in terms of G_. 
Similarly, using i|39[l and (|42|l . we obtain the inverse relation; these read: 



or r U e i(03-0)/2 _ a e -i(,<h-4>)/ 



21 



G + (0!,02,0 3 ) = - Add d0^ 4- G -(■&!, ■&,(/>) ; (47) 

W 7 [0 2 + 2 _ 2 00 2 COS(0 3 - $)] 

or r r#e -i (* 3-< ^/ 2 - ^ne'^ 3 -*)/ 2 ! 4 
G_(0i,0 2 ,0 3 ) = - / d00 / d0 V G+(0i,0,0) . (48) 

W 7 [0 2 + 2 - 2002 COS(0 3 " </>)} 

Hence, if the shear field is a pure E-mode field, these interrelations will be satisfied; note that similar equations are 
valid for the shear 2PCF (Crittenden et al. 2002; Schneider et al. 2002), relating £ + to and vice versa. 

7. Aperture statistics 

In Sect. 5 we defined 3PCFs for the galaxy-mass correlations, and gave practical estimators for measuring them. 
These 3PCFs were then related to the corresponding bispectra through Fourier transform relations in Sect. 6. As was 
demonstrated in Schneider et al. (2004), the relation between the shear 3PCF and the underlying projected mass 
bispectrum is rather complicated and not easy to evaluate. From Eqs. (|37|l through l|40l) we observe that the same is 
true for the relation between the correlation functions considered here and the corresponding bispectra. 

However, we can avoid the use of numerically complicated transformations between the 3PCF and the bispectra 
by considering aperture statistics (Schneider 1996, 1998; van Waerbeke 1998; Crittenden et al. 2002), in the same 
way as has been done for the shear 3PCF (Jarvis et al. 2004; Schneider et al. 2004) and for the measurement of the 
second-order bias factor and galaxy-mass correlation coefficient (Hoekstra et al. 2001, 2002). The main reason for this 
is that the third-order aperture statistics on the one hand provide a very localized measurement of the corresponding 
bispectra, and on the other hand can be readily evaluated from the corresponding 3PCFs. In this section, we apply 
the aperture statistics to the galaxy-mass 3PCFs. The aperture mass is defined as 



(49) 



M ap (0)= J d 2 0lM|tf|) K (tf) = J d 2 0Q e (|tf|) 7 t(tf) .with Qe(#) = ^J d0'0'lM0V W) , 

where Ug is a weight function of zero total weight, i.e., J d00t/e(0) = 0. The scale of the weight function is described 
by 8, and 7t is the tangential shear component as measured with respect to the direction towards the center of the 
aperture. Furthermore, we define the aperture counts 

N{6) = J d 2 U e (\0\) K g (tf) = i J d 2 U {\#\) N{#) , (50) 

where the final equality follows because U is a compensated filter function. We consider here the filter function 
introduced by Crittenden et al. (2002), which has also been used by Jarvis et al. (2004) and Schneider et al. (2004). 
Writing t/g(0) = # -2 it(0/#), the filter function reads 

«(*) - h ( l - e ^ 2/2 ; %) = / d2 - ^ x = £ e "" 2/2 ; Qe{0) - ex P (*) . (51) 



2tt V 2 J ' w/ J " " 2 > -"w 47r6 /4 r y 26> 2 , 

This choice for Ug has the disadvantage that the support of the filter is formally infinite; however, the Gaussian 
factor renders the effective range of support finite, since the filter function becomes extremely small for 3(9. This 
disadvantage over other filter functions that have been employed in cosmic shear studies (Schneider et al. 1998) is 
more than compensated by the convenient mathematical properties, of which we will make extensive use below. 
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We can now relate the third-order correlations between the aperture measures to the corresponding bispectra. 
Consider first 

{M &p {9 1 )M ap {9 2 )N{9 3 )) = {M ap M &p M){9 l ,9^9 3 )=jdH 1 U ei (#i) fd 2 $ 2 Ug 2 {-d 2 ){dH 3 Ug 3 (-d 3 ) (k(?i) K(# 2 )/%(#3 )> ■ 



Replacing the Ki by their Fourier transforms, carrying out the i^-integrations, making use of (|51|l and the definition 
of the bispectra, yields 



d 2 ^ 



(M ap M ap Af)(9 1 ,9 2 ;9 3 ) = J ^-^ J ^^u{txe 1 )u{h6 2 )u{\t l +t 2 \e 3 )b KK& {t x ,t 2 -,-i 1 -i 2 ) . (52) 
Similarly, we obtain for the galaxy-galaxy-mass aperture correlator 

{NMM ap ){9 u 9 2 ;8 3 ) = J ^ J ^H£ 1 e 1 )u(e 2 9 2 )u(\£ 1 +£ 2 \e 3 )b ggK (£ ll £ 2 -,-£ 1 -£ 2 ) . (53) 

Since u is a function that has a very narrow peak, the third-order aperture measures provide very localized information 
on the respective bispectra and are thus ideal for probing the latter. Unless the bispectra have very sharp features, the 
third-order aperture measures contain essentially all information about the bispectra - cf. the corresponding discussion 
for the shear 3PCF and their aperture measures in Schneider et al. (2004). We next show how the third-order aperture 
measures can be calculated directly in terms of the respective correlation functions. For that, it is convenient to 
introduce a complex aperture shear measure, defined as 

M(9):=M ap (9)+iM x (9) = J d 2 tf Q fl (|0|) [ 7t (0) + i 7x (<?)] = / d 2 tf Q fl (|0|) 7 (0; <p) (54) 

where tp is the polar angle of If the shear field is entirely due to the lensing mass distribution, M± vanishes identically 
(Crittenden et al. 2002). A non-zero value for M± would indicate the presence of B-mode shear. In fact M ap (9) vanishes 
identically for B-modes, whereas M±(9) yields zero for a pure E-mode field. Thus, the aperture measures are ideally 
suited to separating E- and B-modes of the shear. 



7.1. Galaxy-mass-mass aperture statistics 

We will now express the third-order aperture measures to the 3PCFs considered in the previous section. First we find 
from the definition l|54|) 

(MMAf) (0i,9 2 ;8 3 ) = J d 2 X x J d 2 X 2 J d 2 Y Q Bl (\X 1 \)Qe 2 (\X 2 \)Ue 3 (| Y\) ( 7 (X 1 ;^ 1 ) 7 (X 2 ;^)« g (F)) , (55) 

where ipi is the polar angle of the vector Xi, i — 1,2. Next we introduce the separation vectors $i = Xi — Y, and let 
(fi denote their polar angle. Then we apply the transformation 

7 (X,;^) = l(Y + 0,; W ) e 2i ^~^ = 7 (Y + 0«; W ) e 2i ^ ^ , (56) 

where we have introduced the notation that a two-dimensional vector X can be written as a complex number X :— 
X\ + iX 2 , which is convenient for expressing phases: If (p is the polar angle of X, then e lip — X/\X\. Inserting the 
foregoing expression into l|55|) and using the definition l|51|) of the filter functions, one obtains 

(MMAf)(9 1 ,9 2 ;9 3 ) = j^J^i J ^ / d ^ ^ G?-(*i,*a) 

x/d 2 y (y* + ^) 22 (V + ^ 2 (2^-|^)ex p r (\ Y ^ ^ l^l 2 ^ l y " 



29j 291 26 *1 



(57) 



The l^-integration can be performed, and the result, when multiplied with the phase factors, only depends on the 
modulus of the "&i and the angle 4> 3 they enclose. Hence, 

/>0O />OC pI'K 

(MMAf)(9 1 ,9 2 ;9 3 )= d^ ^ / d$ 2 d 2 / dfo GL(0i, a , A MMU (Vu ^fo^i, 6 2 ;9 3 ) , (58) 
Jo Jo Jo 
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where the convolution function AmmN is given in the appendix. In the same way, we can evaluate 
(MM*//) (9 1 ,0 2 ;0 3 ) = (47r) 3 g4g4g4 / e2iVl J d ^ e ~ 2 ' 1V2 G+(tfi,tf 2 ) 
x J d 2 Y (f* + 2 (f + tf 2 ) 2 (2^ - |*T) exp 

^oo poo p2tt 

= dtfitfj / dtf 2 tf 2 / d0 3 G+ (^1,^2,03)^^(^,^2, 03 fofc), (60) 

Jo Jo Jo 

where again the convolution function is given in the appendix. By combining (|58J) and Ij59(l . one finds 

(MapMapTV) 2 ; 3 ) = SR [(MMJV) (0 U 2 ; 3 ) + (MM*//) (9 U 9 2 ;9 3 )} /2 
(M±M±/f) (0 1 ,6 2 ; 3 ) = [(MM*//) (0 U 2 ; 3 ) - (MM//) (0 1} 2 ; 3 )} /2 
(M ± M ap /f) (0 1 ,0 2 , 3 ) = 3 [(MM7V) 2 ; 3 ) + (MM*//) (0 X ,0 2 ; 6 3 )] /2 . 

The first of these is the expression that is directly related to the bispectrum fe KKg . The second expression is expected 
to vanish if the shear is caused solely by gravitational lensing. A significant non-zero value of this correlator would 
indicate the presence of B-mode shear which is correlated with the foreground galaxy distribution. The final expression 
is expected to be zero because it is not parity-invariant, as all odd-order statistics in B-mode shear (Schneider 2003); 
hence, a measured non-zero value would indicate that the data violate parity invariance. 

One may consider how the results defined above would change if the modified correlators G± were used in the 
definition of (MM*//) and (MM//). Using the fact that N/N — n g + 1 we find, using Eq. that since U is a 

compensated filter, the additional terms arising from two point shear-shear correlations vanish. Thus we make the 
following useful assertion that when working with the aperture mass statistics, one need not subtract the contribution 
from the shear 2PCF to compute the reduced three point statistics. 



|y+jil 

291 



Y + 2 



291 



291 



(59) 



7.2. Galaxy-galaxy-mass aperture statistics 

Finally, we consider the galaxy-galaxy-mass third-order aperture statistics 

(NNM)(9 X ,9 2 ;9 3 ) = J d 2 X x J d 2 X 2 J d 2 Y Ue 1 (\X 1 \)Ue 2 (\X 2 \)Qe 3 (\Y\) ( Kg (X 1 ) Kg (X 2 ) 7 (y ; $)) , (61) 

where ip is the polar angle of Y . Again introducing the separation vectors = Xi — Y , which have polar angles tpi, 
i = 1,2, using the transformation 

7( y. ^ = 7 fy. fl+^A e i(, 1+ , 2 -2„) = 7 fy. ^_£l\ ^ ^ , (62) 



and inserting the definitions 1|51|) into l|61|) yields 

1 

(^)~9f9l9l 



(////M)(9 1 ,0 2 ;0 3 ) = jjz^mzi / d 2 tf x e 1 ^ / d 2 tf 2 <5(tfr,tf 2 ) 



d 2 Y (29\ - \Y + i?i| 2 ) (29\ - \Y + V 2 \ 2 ) [Y*j exp 



2 



+ \Y + tf 2 \ 2 \Y 



201 291 26»| 



(63) 



dtfitfi/ dd 2 d 2 dfaQWi,Q2,<h)AM M (ti u Q 2 ,fa\0i,ez-,6a), (64) 
Jo Jo 

where the convolution function is also derived in the appendix. The real part of (/T//M) corresponds to the aperture 
3PCF that is directly related to the bispectrum b ggK , whereas the imaginary part vanishes due to parity invariance. 
Finally in this section we note that, as with the galaxy-mass-mass statistics, the quantity (MNM) does not change 
when the definition for Q is used in place of Q in Eq. (|61|l . Second-order galaxy-galaxy lensing effects are therefore not 
important for the three-point aperture mass statistics. 



8. Discussion 

In this paper, we have defined third-order galaxy-mass correlation functions and the corresponding third-order bias 
factor and related them to quantities which can be measured directly through galaxy-galaxy-galaxy lensing. Whereas 
the underlying physical quantities are the bispectra of the mass and galaxy distribution, as well as the cross-bispectra, 
the observables are the shear as measured from the ellipticity of background galaxies in relation to the position of 
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foreground galaxies. We have argued that these third-order correlations have probably been measured already, and 
are almost certainly measurable from data sets existing now; the ongoing and planned wide-field surveys will obtain 
precision measurements of these third-order statistics. The basic observables are the 3PCFs, as their measurement is 
not affected by holes and gaps in the data set which are unavoidable due to the selection of the survey geometry and 
masking of bright stars and galaxies, as well as CCD defects. The 3PCFs are related to the corresponding bispectra 
through a very broad and oscillating filter function (see Schneider et al. 2004 for the case of the shear 3PCF) and 
thus do not directly provide information about the shape of the bispectra. Therefore, we have considered the aperture 
statistics, which on the one hand can be calculated readily from the measured correlation functions, and on the 
other hand are related to the bispectra through a narrow filter function, hence providing localized information on 
the latter. In fact, unless the bispectra vary strongly as a function of their angular wave vectors £i, the third-order 
aperture measures are expected to contain essentially all information about the third-order galaxy-mass correlations. 
Furthermore, the aperture measures allow one to easily identify the presence of B-modes in the shear field, as well as 
unphysical parity-violating contributions, which can only be due to the observing and data reduction steps, unless we 
drop the assumption that the mass distribution in our Universe is parity invariant. 

The convolution functions relating the 3PCFs to the third-order aperture measures have in principle an infinite 
support; in practice, however, the exponential factor implies that the 3PCFs have to be measured only up to a few 
times the corresponding filter scales #j. Therefore, these convolution functions have essentially finite support. 

Nevertheless, in order to calculate the third-order aperture measures, the corresponding correlation functions need 
to be measured to ~ 6# max , where # max is the largest of the three aperture radii. The largest scale on which the 
correlation functions can be measured is determined by the geometry of the survey. For a large survey with contiguous 
area, this scale is limited by the smallness of the signal when entering the linear regime of structure evolution. In 
this case, the aperture measures contain all the information contained in the measured correlation functions. If the 
survey consists of independent patches each of size tp, then the correlation functions can be measured for angular scales 
up to ~ tp, meaning that the aperture measures are limited to # max ip/6. In this case, the measured correlation 
functions contain more information than the aperture measures. In a future work, we plan to study these aspects; 
note, however, that estimates of the 'information content' are tedious, since they involve the covariance of third-order 
statistical measures, and thus depend, in general, on six variables and have to be estimated from sixth-order statistics. 

The correlation function Q measures the same cross-bispectrum as the third-order statistics considered by Menard 
et al. (2003). In that paper, they considered the correlation between pairs of foreground galaxies and background 
QSOs. Due to the magnification of distant QSOs by the matter in which the foreground galaxies are embedded, and 
the steep slope of bright QSO source counts, such a correlation should be measurable with the large QSO sample 
obtained by the Sloan Digital Sky Survey. Instead of magnification, we employ the shear around galaxy pairs. Which 
of the two methods is better able to measure the galaxy-mass bispectrum depends on the available observational data. 
For deep wide-field images, the shear method investigated here will be more efficient, given the sparseness of bright 
QSOs on the sky. 

As is true for the second-order galaxy-mass correlations, as measured in GGL, the physical interpretation of the 
third-order correlations is not straightforward, but needs to be done in the frame of a model. On very large scales, we 
might expect that 63 essentially becomes a constant, but that will be difficult to measure, as on these large scales, the 
density field is expected to quickly approach a Gaussian, and thus the bispectrum should be very small. 

A useful analytic description of clustering statistics on all cosmological scales is given by the halo model (see Cooray 
& Sheth 2002 for a review). Within this framework it will be possible to relate the higher-order cross-correlation 
functions to the quantities that specify the Halo Occupation Distribution. As pointed out by Bcrlind & Weinberg 
(2002) the HOD essentially contains all of the information about the statistics of galaxy clustering that theories of 
galaxy formation are able to provide; it is a complete description of the bias. Future empirical determinations of the 
HOD using observations of galaxy and mass clustering will enable powerful constraints to be placed on the development 
of models for galaxy formation. In a forthcoming paper we shall explore these concepts in detail, building the physical 
interpretation of our third-order galaxy-mass correlations using the halo model and the HOD. 



Acknowledgements. We thank Martin Kilbinger for useful comments on the manuscript. This work was supported by the 
German Ministry for Science and Education (BMBF) through the DLR under the project 50 OR 0106 and by the Deutsche 
Forschungsgemeinschaft under the project SCHN 342/3-1. 



Peter Schneider and Peter Watts: Galaxy-galaxy-galaxy lensing 

Appendix A: The convolution functions 



13 



In this appendix, we derive the convolution functions which relate the third-order aperture measures to the corre- 
sponding 3PCFs. We start by rewriting the exponent in l|57|l. (|59|l and 1)63(1 by making a translation Y = y — c, where 
c is chosen so as to remove linear terms in y, 



29\ 



\y + # 2 \ 2 

291 



|!T_|yK 



9fl2/)2/)2 



a 2 / ■&! i?: 
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(A.l) 



h?i| 2 

29\ 



H ®2 «2 /0f 2^2 COS 03 ^ 



a 2 



2B\ 29\ 



9f 



+ 



+ 



(A.2) 



Next, we consider the prefactor of the exponential in (|57|) . and change the integration variable from Y to y. For this, 
wc first define 



for i = 1, 2, or more explicitly. 

a 2 /^i i? 2 e" i</>3 



c e 



.91 = 01 



; 32 = ^2 - ^ 7JJ + -^5— 



(A.3) 



(A.4) 



The prefactor of the exponential, when multiplied with the phase factors, then becomes 

KimM - (y*e ilfil + 9if (y*e iV2 + hf (2^ 3 2 - \V\ 2 - |c| 2 + cy* + c*y) 

= g\g\ (291 - \v? - l<f) + 2|y| 2 3iff 2 [h <p\ - 9i) + 9i 0*a - &)] + C = F UMN + C, 



(A.5) 



where C denotes additional terms which, however, become zero when the y-integration is carried out, since the 
exponent depends only on \y\ 2 . The y-integration can now be carried out, 



m ma/" c 12/1 /a2 = 7ra^!ff 2 



The result depends only on d\, ?9 2 and the angle <j) 3 between i?i and i9 2 , so that an additional angular integration can 
be carried out in Eq. 1(57(1 . We thus arrive at 1(58(1 . with 



Amaiat 



9ih e~ b ° 
72ire 8 



7# + -£ ~ — 5132 + 2 (ff 2 i?i + 5ii? 2 - 25i3 2 ) 
t>i t/ 2 a 2 



e 4 



g 2 /) 2 i fl 2 fl 2 i Q2r\2 



(A.6) 



Similar expansions of the prefactors of the exponentials in Eqs. 1(59(1 and ((63|) can also be determined using the 
same definitions as above. Performing the relevant integrals, we find that the convolution factors for equations 160(1 
and ((64(1 are then given by 



-b 



Amm*N 
and 



727t8 8 
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( (Si - #i) (32 - #2) [ — FiF 2 - (Fi + F 2 ) + 2a 2 + gj&e -1 * 3 
I L a 2 
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(A.7) 



(A. 



where F, = 29j 
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